Unified Theory of Interspecific Allometric Scaling 
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A general simple theory for the interspecific allometric scaling is developed in the d+l-dimensional 
space (d biological lengths and a physiological time) of metabolic states of organisms. It is assumed 
that natural selection shaped the metabolic states in such a way that the mass and energy d + 1- 
densities are size-invariant quantities (independent of body mass). The different metabolic states 
(basal and maximum) are described by considering that the biological lengths and the physiological 
time are related by different transport processes of energy and mass. In the basal metabolism, 
transportation occurs by ballistic and diffusion processes. In d = 3, the 3/4 law occurs if the ballistic 
movement is the dominant process, while the 2/3 law appears when both transport processes are 
equivalent. Accelerated movement during the biological time is related to the maximum aerobic 
sustained metabolism, which is characterized by the scaling exponent 2d/(2d -I- f) (6/7 in d = 3). 
The results are in good agreement with empirical data and a verifiable empirical prediction about 
the aorta blood velocity in maximum metabolic rate conditions is made. 

PACS numbers: 87.10.-(-e, 87.23.-n 



Metabolic rate B and body mass M are connected by 
the relation B — aM^ , where h is the allometric exponent 
and a is a constant. For several decades it was accepted 
that the basal metabolic rate (BMR) among almost all 
organisms [J [2i @| (interspecific scaling) was character- 
ized by & = 3/4 (Kleiber's law Q). A few years ago, 
theoretical explanations of the ubiquity of the 3/4-law 
based on the resource distribution network common to 
all organisms 0] and on network geometry optimization 
d, 0] were proposed. Kleiber's law has however been 
questioned recently. On the observational side, it is not 
clear whether the value of & is 3/4, 2/3 or even variable in 
both interspecific and intraspecific (same species) scaling 
d, H, On the theoretical side, there are several de- 
bates about the validity of these models ^, TT*! undermin- 
ing the rational basis for the scaling law. A related open 
question is why maximum aerobic sustained metabolic 
rate (MMR) of endothermic animals scales with an ex- 
ponent larger than that of BMR d, [13, [ill, [13, [11 El ■ 
In this work we show how, on the basis of a few quite 
general hypotheses, all the aspects of metabolic scaling 
mentioned above may be accounted for. The central as- 
sumption of our approach is that, if one characterizes 
the metabolic state of an organism as a point in a space 
comprised of d biological lengths and a biological time 
T, then the {d -\- l)-dimensional mass and energy den- 
sities are size-invariant under natural selection. Using 
this idea we may, given the dominant mass and energy 
transport processes, obtain dynamical relations involving 
lengths and time. The latter are characterized by invari- 
ant quantities (diffusion coefficient, velocity etc.). Each 
transport mechanism is related to a different metabolic 
state, with its own values for allometric exponents, in 
agreement with data. 

The key questions in interspecific allometry are: (1) Is 



there a universal interspecific BMR exponent 6? (2) Is 
there a model that can describe both BMR and MMR 
scaling? (3) Which quantities have their scaling deter- 
mined by the BMR, and which by the MMR? 

We present a simple and general theory that answers, 
at least partially, these questions and explain some as- 
pects of the allometric scaling of organisms. As all biolog- 
ical processes depend on characteristic times, it is natural 
to include a characteristic time, as well as various charac- 
teristic lengths, when specifying the metabolic state of an 
organism [1^, [l^| . We therefore associate the metabolic 
state of an organism with a point in a d + 1-dimensional 
space with d biological lengths Li, L2 . . . , Ld and a bi- 
ological time r. Here d is the number of spatial dimen- 
sions; although we usually have d = 3, we will work in 
general dimension d. Examples of biological lengths L 
and times r are the total aorta length in mammals, the 
length of capillaries, the mean distance from cell surface 
to mitochondria in unicellular organisms, the duration 
of one heartbeat, the capillary blood transit time or the 
turnover time for glucose P, [3, [3, [l3| ■ 

We now look for simple, general relations constrain- 
ing the distribution of points in the space of metabolic 
states. It is usual for an animal to make several transi- 
tions between states A^asai and Amax, through a series 
of complex biochemical processes. We do not consider 
these transitions. Rather we are interested in describing 
how natural selection has shaped the state A^asai and 
other ones belonging to the space of metabolic states 
of all organisms (or a group of them). An organism 
is characterized by the fundamental quantities of mass 
and available energy for the metabolic processes. We 
therefore identify the mass M and available energy E 
as the fundamental variables characterizing an organism. 
Since we have a c? + 1-dimensional space we shall use the 
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mass density pd+i{Li, . . . , Ld, t) (mass per unit volume 
and unit time) and the energy density ad+i{Li, . . . , L^, t) 
(available energy per unit volume and unit time). We as- 
sume that during evolution, natural selection enforces the 
constraints of size-invariant (independent of body mass) 
Pd+i and (Jd+i (our first and second hypotheses, respec- 
tively). Our third hypothesis is that the scaling of the 
metabolic states is determined by the dominant dynam- 
ical transport processes of nutrients (mass and energy), 
which are characterized by size-invariant quantities (dif- 
fusion coefficient, average velocity, etc.). Note that the 
first and second hypotheses furnish two relations valid 
for all metabolic regimes. Different metabolic scalings 
will appear because there are different ways to transport 
nutrients. 

From the second hypothesis we obtain that E — 
fd+i T^d, where rV^ is the d + 1-volume and Vd = 
L1L2 ■ ■ ■ Ld- Since power is defined as P = dE/dt, energy 
can be written in terms of B, the power averaged over the 
time scale t, as E — Bt. We identify B as the metabolic 
rate. Therefore from the first and second hypotheses we 
obtain that 

M = pd+i rVd, and B = (Jd+i Vd ■ (1) 

Let us briefly present a qualitative argument about an 
optimal delivery that supports our postulates. Consider 
two organisms with the same body mass M belonging 
to the same group. The organism with the larger bi- 
ological volume Vd has the larger nutrient distribution 
network and more fuel and oxygen are arriving to the 
cells in a unit time. Therefore its cells must have a fast 
metabolism in order to consume the fuel. Of course, a 
fast metabolism is related to a small biological time. On 
the other hand, a small Vd implies in a larger biological 
time T. These arguments suggest that the product rVd 
is constant for these organisms. 

An immediate consequence of our hypotheses is that 
the power per mass (specific metabolic rate), namely 

B_ _ (7d+l 1 
M pd+l T 

is inversely proportional to the metabolic time. Animals 
with a small r, such as small mammals, require larger 
power by unit mass than ones with a large r (large mam- 
mals) because their cells have a large mitochondrial den- 
sity pj]. 

Obviously, r cannot be neither 00; there must be 
a minimal and a maximal metabolic times Tmin and 
Tmax- For auimals, a lower bound for Tmin can be found 
from the observation that the biological volume of an 
organism Vd cannot be larger than its spatial volume 
V, since the biological lengths Li characterize the or- 
ganism's anatomy on the scale of the body, or some 
organ or cellular structure. For compact animal bod- 
ies we have that M = pV, where p is the usual d- 
dimensional mass density, which is approximately con- 
stant. Using pV = M = pd+i TVd in > Vd, we obtain 



T > Tmin = p/pd+1- Siuce B/M cannot be zero it must 
exist an size- invariant {B /M)min related to the mini- 
mum power per mass to keep the organisms just alive, 
the so-called tissue maintenance specific metabolic rate 
[13]. This implies that T„iax = ad+i/{pd+i{B/M)jnin)- 
Note if T is size-invariant, such as Tmax, we have isomet- 
ric scaling of the metabolic rate (5 = 1). 

The relation M oc rVd, derived from the first hypoth- 
esis, is a generalization of the result of Banavar et al. 
0, namely Vnet oc L'^'^^, where Vnet is the total vol- 
ume of an efficient distributive network. Using that the 
blood volume Vnet is proportional to mass, they obtained 
M (X L''"'"^, a basic relation to deduce the 3/4-law. The 
two relations are equal when t (x L, a condition valid for 
the BMR. Moreover, the result M cx L''+^ is also cru- 
cial to obtain the BMR exponent in the model of West, 
Brown and Enquist (WBE) This relation is also a 
generalization of the equation F oc {Lp/u)B of Banavar 
et al. Q, where Lp is the physical length of the system 
and u is the characteristic length scale. If we rewrite this 
equation as M oc ptVd, where p is the tissue density and 
t is the physiological time related with the rate of energy 
use per unit volume, it becomes similar to Eq. ([T]). 

Although, from our third assumption, we need some 
dynamical size-invariant quantities, like the blood fiow 
speed velocity vo in the aorta or in the capillaries, the 
length Ic and the radius Vc of capillaries are not neces- 
sarily invariants. 

Let us first study the case of transport via diffusion. 
We have only one metabolic length scale (Li oc L2 oc 
. . . oc Ld oc L), so the biological volume is given by 
Vd oc L'^. Since diffusion over short distances is fast, 
it is possible that the metabolic rate of very small organ- 
isms is governed by this process. It is well known that 
L = Dqt^/'^, where Dq is the size- invariant diffusion co- 
efficient. Since r = {L/D^Y, we obtain from Eq. ^ 
that M oc L'^^'^. This relation furnishes how L depends 
on M and we can use again Eq. ([1]) to obtain that 

L oc , T (x , 

B oc Af^fd . 

In d 3, the metabolic exponent is b = 3/5. 

For larger organisms diffusion is inadequate. Transport 
by convection is then utilized on large length scales. In 
mammals, for example, we find the cardiovascular system 
that transports blood to the capillaries, where the cells 
are fed by diffusion. Since blood circulates in an ballistic 
regime, we consider that the BMR is basically driven by 
ballistic transport, namely L = vqt, where the velocity vq 
is size-invariant. Now we must specify how the different 
metabolic steps are related. We call BMR-1 the scenario 
of a single metabolic relevant length Li (x L = vqt and 
a single time r, both related to the ballistic transport. 
The other lengths, related to other metabolic steps, have 
evolved to meet it, namely L2 (x . . . oa Ld (x L. Using 
that Vd oc L'^ in Eq. ([T]), we write 

L oc r oc M rfTT 
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B oc . 

For d = 3 we find the 3/4-law, namely t oc L oc Af^/* 
and B oc M^^^. These results are the same as those of 
WBE [i] and Banavar et al. [i,[3]. 

In the BMR-2 scenario, we have different relevant 
lengtlis and times related to the metabolic processes. 
However, due to the concept of symmorphosis [ig, which 
states that all metabolic steps have co-evolved in order 
that no step is more limiting than another, we will end 
up with a single time r and d — 1 rescaled lengths. In a 
"cylindrical" symmetry we have Li oc L — voti (ballis- 

1 /2 

tic term) and d — 1 lengths proportional to i? = -00^2 
(diffusion). Both vq and Dq are size- invariant. From 
the symmorphosis principle [ti — t2 = t), it follows 
that R = {Do/vq^^)L^/'^. The biological volume is 
Vd oc R'^-^L. From Eq. (P) we obtain that 

L oc T oc , R oc ^ 

l + d 

B oc M 3+<i . 

Then in d = 3 the BMR-2 scenario yields the 2/3 law, 
without mention of the area/volume ratio. We obtained 
both 3/4 and 2/3 laws from the same transport processes: 
convection and diffusion. If convection is the dominant 
limiting process we have the 3/4 law; if the two processes 
are equivalent we obtain the 2/3 law. 

The circulatory networks of endothermic animals are 
dynamical ones which are adjusted according to the 
metabolic state. The transition from resting to maxi- 
mum activity can be described very briefly as follows: 
(a) the heart increases its rate and output; (b) the mean 
arterial pressure and peripheral extramuscular resistance 
increase; (c) arterial blood volume increases due to con- 
striction of the veins; (d) extramuscular flow remains 
essentially constant, somewhat reduced in some organs 
but increased in others; and (e) total flow and mus- 
cular flow increase, with all muscular capillaries acti- 
vated. The items (a), (c) and (e) suggest that we have a 
"forced movement" during the characteristic time r. This 
means that the typical constant velocity can be written as 
V = aoT, where oq is a size-invariant acceleration. Con- 
sequently the MMR is driven by an inertial movement 
accelerated during time r, implying that L = vt = qqt^ . 

If inertial transport is the only relevant process (MMR- 
1 scenario), it follows that Li oc L2 ■ ■ ■ oc Ld oc L = aqt^. 
Since Vd oc L'^ and r oc L^/^, we obtain from Eq. ([T|) the 
metabolic relations: 

L oc , T oc M^arr ^ 

B oc M53TT . 

For d = 3 we have that L oc M^/'^, r oc M^/'^ and B oc 
M^/^. This results agree with the ones obtained trough 
a generalization of WBE ideas to MMR scenario ■ 

In the MMR-2 scenario, diffusion and inertial move- 
ment are equally relevant. We again choose a "cylindri- 
cal" symmetry so that have that Li oc L = apT^ while 



TABLE I: Allometric exponent y describing the dependence of 
a variable Y on body mass M (F ~ M^). Under parenthesis 
is the error in the last significative of the observed quantities. 
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the remaining d—1 lengths are of order Dqt^^^. It follows 
from Eq. ^ that 

4 2 
L oc M ti+s , r oc M'i+s , 

B oc , R oc . 

When d = 3, we obtain that L oc M^/^, r oc M^/^, 
R oc A/i/s and B oc M^/^. 

The transportation network can be characterized by 
"aorta" and "capillaries". Note that the aorta La and 
capillary Ic lengths are both proportional to L. Since 
the nutrient fluid is conserved, the volume rate of flow is 
given by 

Q = nRlva = NcTrrlvc , 

where Ra and Va are the aorta radius and fluid velocity 
and Nc, Vc and Vc are capillary number, radius and fluid 
velocity, respectively. It is natural to write that Q oc B. 
In the basal regime, Va and Vc are size-invariant. Then 
we obtain that Ra oc B^^^ and NctI oc B. Making the 
extra assumption that Tc is invariant, it follows that the 
capillary density — N^/M is Pc oc B/M oc l/r. For 
d = 3 we have the following results: (a) BMR-1 - Ra oc 
Af3/8 and Pc oc M'^/^; (b) BMR-2 - Ra oc Af^/^ and 
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Pc oc M~^/^. In the maximum regime, we have that 
Va — o,qt, implying that now Va depends on the mass. 
Note that this suggests new empirical studies. Then, 
we obtain that i?^ oc B/t. Now Vc is not necessarily 
invariant, as in the basal case. Since pcT^Vc oc 1/t we 
obtain that pc cc 1/t if we make the extra assumptions 
that both Vc and are independent of body mass. The 
results are: (c) MMR-1 - Ra oc M^^^^ and pc oc M'^/''; 
(d) MMR-2 - Ra oc Af and pc oc Af-i/*. 

Let us compare our predictions for d = 3 with exper- 
imental data (see Table The values predicted in the 
MMR-2 context are far from the experimental ones. On 
the other hand, the MMR-1 scenario describes very well 
the MMR data. The exponent 6 = 6/7 « 0.86, larger 
than the basal value, is in very good agreement with 
data. Muscular capillary density of mammals is linked to 
MMR, instead of BMR, because only during exercise are 
all the muscular capillaries perfused. The capillary den- 
sity scales as Pc — N^/M oc M~^^'^ , in good agreement 
with the ave rag e experimental value for various regions 
of muscle Since r oc M^/'', frequencies must scale 

as F oc oc Mf with / ^ -1/7 w -0.14. This value, 
smaller than the basal one, is also in good agreement 
with data for heart and respiration rates in strenuous 
exercise. The results for the capillary invariant radius 
Tc and Ic oc M^/^, agree roughly with the theoretical- 
empirical estimation of Dawson {Ic oc M^-^^ and 
re oc MO-08). 

Since the length of aorta cannot change from basal 
to maximum metabolic regimes, it should scale as the 
prediction of MMR scahng. The predicted exponent 0.29 
agrees well with data. The aorta radius could in principle 
follow the two scalings because of the elasticity of the 



aorta and the dynamical body adaptations of mammals 
in the transition BMR - MMR. The experimental value 
0.36 has however a better agreement with the MMR-1 
value 5/14. 

Consider now the predictions of BMR scaling. Re- 
cently the empirical values of the BMR exponent of mam- 
mals were reanalyzed using diverse procedures [1, 
that furnished values in the interval between 62 = 2/3 
and 61 = 3/4, which are the predicted values of BMR- 
2 and BMR-1, respectively. Heart and respiration rates 
are close to the BMR-1 value —1/4 and other biological 
variables have values close to multiples of 1/4. On 
the other hand, empirical data near multiples of 1 /3 are 
also reported [lOl]. Therefore, the two scenarios are pos- 
sible. This last possibility explains why b is greater in 
large versus small mammals data: diffusion and ballistic 
transports can be equally important in small organisms 
(BMR-2) but not in large ones, where ballistic transport 
is crucial (BMR-1). Finally, let us emphasize that we 
make a verifiable empirical prediction: the aorta blood 
velocity (wq), which is scaling-invariant in BMR condi- 
tions, grows with mass in the exercise-induced MMR 
condition {va oc r). The related exponent, which is pre- 
dict to have the value —0.14 in the MMR-1 scenario, was 
never measured. Its empirical determination can be an 
experimental test of the importance of the transportation 
processes for the allometric scaling of metabolism. 
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